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ABSTRACT 

Empirical models for emission from stars and the ISM are subtracted from the zodiacal- 
light-subtracted DIRBE 3.5 /un emission. Because the models are contaminated by unknown 
levels of the CIB at other near-IR wavelength, the residual is not simply the 3.5 /im Cosmic IR 
Background, but a linear combination of the background levels at several wavelengths. In spite 
of this, the residual can be used to place limits on the near-IR CIB intensity if its spectral shape 
is assumed. Additionally, the residual level is shown to be more nearly isotropic than previous 
estimates over a much larger fraction of the sky. An excellent correlation of near-IR and far-IR 
ISM emission provides evidence of the high accuracy of the brighter stellar emission model. The 
possibility that any residual emission is zodiacal in nature is discussed. 


1. Introduction 

TBD: The introduction will describe the Cosmic Infrared Background (CIB) and its significance. It will 
describe prior efforts at its detection, and point out successes and limitations, setting the stage for the new 
modeling presented here. Hauser et al. (1998), Arendt k Dwek (1998), Gorjian, Wright, k Chary (2000), 
Matsumoto, et al. (2000), Kashlinskv. et al. (1996a, 1996b), Kashlinsky k Odenwald (2000). 


2. The Models 

At the far-IR wavelengths Hauser et al. (1998) modeled the ISM emission at 140 and 240 pm by 
assuming that it was linearly correlated with the 100 pm emission that was used as a template. This worked 
moderately well at high latitudes, but left obvious defects especially at low latitudes. The reason is that 
that linear scaling of the 100 pm ISM template by a single constant coefficient assumes that all of the ISM 
emission has the same spectrum. Thus regions warmer or colder than average show up as positive or negative 
residuals. An improved model was found by noting that the ISM emission followed a well-defined trend in 
color-color plots. The narrow trend implied that the 240 pm ISM emission could be modeled as a linear 
combination of the 100 and 140 pm emission. The differences between the two ISM templates would account 
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for variations in the spectrum of the ISM even when each template is scaled by a single coefficient. Indeed, 
the isotropy of the residual emission after subtraction of the two-component model was improved enough 
over a wade range of latitudes to justify the claims for detection of an isotropic CIB. 

At near-IR wavelengths, the Dwek k Arendt (1998) approach for modeling the Galactic stellar emission 
was that of a one-component model. The 2.2 //in emission w r as scaled by a single coefficient in order to model 
the stellar emission at each of the other near-IR bands. This empirical model proved to be more accurate 
than the statistical model used by Hauser et al.(1998), both on the small scales where the statistical model 
could not match the actual placement of individual stars, and on the large scales, where the statistical model 
did not properly match the Galactic intensity gradient. However, like the one-component ISM model, the 
Dwek k Arendt model cannot properly subtract the stellar emission in regions where the mean spectrum 
changes, caused either by extinction or by changes in the intrinsic colors of the stellar populations. Thus, 
we seek an improved model of the stellar emission by using a two-component model. 

Our basic approach is to model the stellar emission at 3.5 /im as a linear combination of the emission 
at 2.2 and 4.9 //m: 

/.,»r(3.5) = A/ Star (2.2) + B/ S , ar (4.9). (1) 

We refer to this model of the 3.5 /mi emission as “Model 1”. The model appears to work well at high 
latitudes, but runs into problems when we try to apply it at low latitudes where extinction affects the stellar 
colors. Thus, to correct for extinction effects, with the intent that the coefficients .4 and B will then model 
intrinsic color variations in the skv brightness, we use the more complicated expression (“Model 2’ ): 

W(3.5) = .4/.i.r(2.2)e- (T »- < -™ ) + BI star (4.9)e^~ T ^ ] . (2) 


The extinction correction factors are estimated using the fact that 


/star (I -25) e -(n 25-T2.2) (3) 

/star (2-2) 


This ratio of the two shortest wavelength bands is used because extinction will be strongest, and ISM emission 
appears to be negligible in these bands. Thus, the actual model for the stellar emission we apply is given by 

~ / star (1.25) l J (4 , 

I, tar (2.2) J 


/star (3.5) = ,4/, tar (2.2) 


' [staA 1-25) • 

[ W( 2-2) . 


+ B/, (nr (4.9) 


where a = (r 3 . 5 - r 2 . 2 )/(r 1 . 25 - r 2 . 2 ) and 3 = (r 35 - u. 9 )/[t im - n.t). The numerical values of a = -0.317 
and 0 = 0.206 are calculated using the Rieke & Lebofsky (1985) reddening law, which has been shown to be 
a good representation of the extinction observed by DIRBE (Arendt et al. 1994). 

ISM emission is present at 3.5 ^tn, though very difficult to see without subtraction of the stellar emission. 
Our model of the ISM emission will be the simple one-component model used by Hauser et al. (1998): 


f/SM(3- 5) = C7 /sa/(100) 


( 5 ) 


where the 100 /xm ISM intensity is the observed emission after subtraction of a 20 nW m 2 sr 1 estimate of 
the CIB at that wavelength (Hauser et al. 1998). 

The final step is to combine these models with a constant term, D (potentially the isotropic CIB), to 
model the observed 3.5 /xm emission as 


U,(3.5) = AI ob ,C2.2) + BI ob M-9) + CJ, s .w(100) + D 


( 6 ) 
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for Model 1, and 


Jots(3.5) = AI 0 b 9 {2.2) 


r \ 1 J 


[W 1 - 25)1 

f* 

+ BUA 4 . 9 ) 

7 ^ ( 1 . 25 )' 

T, 6 s( 2 . 2 ) j 


lobs ( 2 - 2 ) 


C I ISM (100) + D 


for Model 2. The coefficients .4, B, C. and D ran now be determined with a simple least-squares fit. In 
deriving the solution we excluded the same bright stars and low latitude emission that were excluded at 3.5 
Itm in the Hauser et al. (1998) analysis. For Model 1, we also found it necessary to exclude regions at low 
Galactic latitude (|6| < 10°), in order to prevent the model from attempting to match low latitude color 
variations at at the expense of a good fit at high latitudes. The derived values of the model coefficients are 
listed in Table 1. In Figure 1. we show the residual map that results from subtraction of the stellar and 
ISM components of each model from the 3.5 micron emission. This residuals have a mean level of D (in the 
unblanked regions), and are depicted on a comparable intensity range as the Hauser et al. (1998) and Dwek 
& Arendt (1998) residuals that are also shown in Figure 1. 


Note that w r e actually use the observed 1.25, 2.2, and 4.9 //m intensities in these models rather than 
only the purely stellar components, since the levels of the CIB at these wavelengths are to be determined. 
The consequences of this are important, and are discussed in the next section. 


3. Interpretation of the Residual 

Ideally, if we could have used templates of stellar emission that were free of CIB emission at 1.25, 2.2, and 
4.9 fj , m and free of ISM emission at 4.9 gm. then the residual emission should represent the CIB at 3.5 fi m. 
Instead, these unsubtracted backgrounds propagate through the model and affect the level of the residual. 
For Model 1, the effects can be traced by rewriting Eq. 6, explicitly identifying all of the components. 

D = [/s<«r(3.5) + //sm(3-5) 4- 7c/b(3.5)] — 

A[I stnr {l.2b) + //sm(T 25) 4- /o/b(T 25)] - 
B[I sta r( 3.5) + //sm (3-5) + /c/a(3.5)] - 

Cl ism ( 100) W 

Next, if we assume that the stellar and ISM models are good {i.e. I s t ar (3.5) — A/star (2-2) 4- B I star (4.9) 
and W(3.5) = B//sa/(4.9) 4- CIism (100)}, then the constant term D can be expressed as a function of 
the CIB intensity at several wavelengths: 

D = /c/s(3.5) - AIcib{ 2-2) - BIcibU- 9). (9) 

When we apply the coefficients derived for Model 1, this expression becomes: 

D = Icib{ 3.5) - 0. 445/c;/ B (2.2) - 0.220/ C /s(4.9) = -0.00106 (10) 

where all intensities are in MJy sr' 1 . or 

D = vIcib( 3.5) — 0.279 vIcib{ 2.2) — 0.308 ^/c/b(4-9) = —0.91 (11) 

with the vl v in nW m -2 sr“T 

For Model 2. the extinction correction makes the equivalent expressions more complex, but with some 
algebraic reshuffling and keeping only first-order terms, we find 

D - Icib{ 3.5) - [AarfT 1 4- £3rf 2 r 4 i]/c/B(1.25) - 

[A( 1 - a)r“ 2 4- B/3rf 2 r 42 ]/c/B(2*2) - Br^ 2 I C iB{ 4>9) 


( 12 ) 
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where n 2 = / s(ar (1.25)// S („ r (2-2). etc. Inserting the numerical values for r 12 = 1-2. r 41 - 0.25. r 42 - 0.3 
and for the other terms, we finally arrive at 

D = Icib{3-5) + 0.103/r/fl( 1-25) — 0.564 Iojb(2.2) — 0.259/c/b(49) = -0.00234 (13) 

where all intensities are in MJv sr _1 . or 

D = vIcib{ 3.5) + 0.0386 i//r/ fl (1.25) - 0.355 vIcib{ 2-2) - 0.363 x//c/b(4.9) = -2.01 (14) 

with the vl u in nW m " 2 sr 1 . 

The above expressions show that for both Model 1 or Model 2, the unsubtraeted backgrounds in the 
2.2 and 4.9 /xm components of the models reduce the value of the residual D from the level of the 3.5 /xm 
C1B. Conversely, the 1.25 ^rn component adds a small amount to the residual level in Model 2. If the 100 
^m CIB had not been subtracted from the template of the ISM emission it would also reduce the observed 
residual levels by -0.0387 t'/c/fiU 00) nW m -2 sr -1 in Model 1, or -0.0347 b/c/bO 00) nW m -2 sr -1 in 
Model 2. 

As discussed by Hauser et al. (1998), the systematic uncertainties dominate the random statistical 
uncertainties in obtaining the mean 3.5 /xm residual. Systematic uncertainties include the detector offset 
and zodiacal light uncertainties, and the uncertainty in the stellar and ISM model developed here. The 
zodiacal light uncertainties used by Hauser et ah (1998) are equally appropriate here, and are reproduced 
in Table 2. Detector offset uncertainties are negligible in comparison, and are therefore not listed. The 1 o 
uncertainties in the new stellar and ISM models are estimated by using the full variation in the mean residual 
intensity of the various high latitude regions listed in Table 3 (see below). As in Hauser et al. (1998), we 
take the quadrature sum of the zodiacal light and combined stellar and ISM uncertainties to arrive at the 
total systematic uncertainties, listed for each model in Table 2. The scale factors that are applied to the 
zodiacal light uncertainties are the same as those used to scale each wavelength in the construction of the 
stellar and ISM models. 


4. Isotropy of the Residual 

While the residual emission derived here after subtraction of the 3.5 /im stellar and ISM models is not 
the CIB directly, it should still be isotropic if the models are accurate. Therefore, following the example of 
Hauser et al. (1998), we have examined several different means of assessing isotropy. In all cases, the key 
issue is whether or not the data show only the amount of variation expected wdthin the uncertainties. 


4.1. Mean Patch Brightness 

The simplest test of isotropy is comparison of the mean brightness in selected patches. Table 3 lists the 
brightnesses derived at several of the regions examined by Hauser et al. (1998) and Dwek & Arendt (1998). 
The present variation in brightness among these patches is smaller than in either of the two previous studies, 
and is smaller than the expected systematic uncertainty (even if only the zodiacal light component of the 
uncertainty is considered). The patch intensities for Model 1 are not quite as uniform as those for Model 2. 
but in both cases the variations smaller than the expected systematic uncertainties. 



4.2. Brightness Distributions 


The next isotropy test is the examination of the intensity histograms of the residual emission. For a 
truly isotropic residual, the shape and width of the histograms will reflect only the random noise uncertain- 
ties. Histograms constructed for the north and south portions of the high-latitude HQA and HQB regions 
examined by Hauser et al. (1998) are shown in Figures 2 and 3. The results of fitting Gaussian distributions 
to the histograms are shown in the figures and listed in Table 4. The new results exhibit smaller dispersion 
than the previous results, particularly over the relatively large HQA regions. North-south asymmetries in 
the means of the distributions are also reduced in the new results. Kolmogorov-Srnirnov (K-S) tests applied 
to the distributions, indicate that the HQBN residuals are essentially Gaussian for the new results with 
Model 1 and especially Model 2. while distributions in the HQBS (and for the Dwek k and Arendt (1998) 
HQB regions) are probably only slightly non-Gaussian. The K-S probabilities of the distributions differing 
from Gaussian are listed in Table 4. A probability of 1.00 is definitely non-Gaussian. a probability of 0.0 
indicates a perfectly Gaussian distribution. 


4.3. Systematic Spatial Variations 

The tests above are necessary conditions for isotropy, but they are insensitive to any anisotropic structure 
that does not alter the mean level or introduce some skewness or tails to the brightness distributions. A 
simple means of looking for likely spatial structure is to check for gradients in the residual emission. Figure 
4 shows the gradients in the residual emission as a function of the cosecant of Galactic latitude. Data for the 
north and south hemispheres are shown separately. By this test, the Dwek & Arendt (1998) model removed 
Galactic emission more effectively than the models used by Hauser et al. (1998). The present models are 
seen to make further improvements, most dramatically at lower latitudes. The residual emission of Model 
2 exhibits smaller gradients than that of Model 1 over this range of latitudes. Table 5 lists the gradients 
for the HQA and HQB regions, showing no significant Galactic gradient in the HQB region for the residual 
emission derived from Model 2. The derived gradients in the HQA region are even smaller, but statistically 
significant as the region is about ten times larger. 


4.4. Two-Point Correlation Functions 

The most rigorous test of spatial isotropy that was used by Hauser et al. (1998) is the two-point 
correlation function. A truly isotropic residual should show no structure above the expected random noise 
on any spatial scales. The two-point correlation functions for the 3.5 residual emission from the present 
and previous results are shown in Figure 5. Note that because the Hauser et al. (1998) results relied on a 
statistical model for the stellar emission, the mean level but not the detailed structure of the faint stars was 
removed from the data. Thus, the Hauser et al. results have a much higher “random noise level than either 
the Dwek et al. (1998) or the present results. The variance of the residual emission in the present results is 
slightly smaller than for the Dwek k Arendt results. The two-point correlation function over the HQB region 
for Model 1 appears to be very nearly random by this test. The Model 2 results are not quite as isotropic as 
the Model 1 results, but a distinctly better than the Dwek k Arendt results. Detailed examination of the 
images (Figure 1) suggests that the 2 point correlation functions in the HQB regions are sensitive to small 
systematic errors in the structure of the zodiacal light model. With the two-component models used here, 
errors at different wavelengths partially cancel one another. This cancellation turns out to be slightl\ more 
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effective for the parameters of Model 1 than those of Model 2. 


5. The Interstellar Medium 

A good model of the Galactic stellar emission is necessary tool in the study of the relatively faint 
Galactic ISM at 3.5 and 4.9 /xm. Freudenrcich (1996) shows that emission of the ISM can be be traced to 
high latitudes using the 3.5/2. 2 /mi colors. Arendt et al. (1998) correlated the variations in the reddening- 
free near-IR colors with 100 /ini ISM emission to derive the 3.5/100 /mi and the 4.9/100 /xm colors of the 
ISM. However, the derivation of these colors was restricted to very low galactic latitudes. Now, with the 
stellar emission models developed here, we can subtract the stellar emission directly, producing a map of the 
3.5 /x m ISM emission instead of only a near-IR color. The 3.5 /xm ISM emission derived from subtracting 
the stellar emission of Model 2 is shown in Figure 6. Comparison with the 100 micron emission reveals an 
excellent match of features down to the noise level of the 3.5 /xm map. The linear correlation of the 100 and 
3.5 /im emission is plotted in Figure 7 for regions |6| < 20°. The linear trend plotted is that derived from 
Model 2. The inset of shows that the linear trend fits all the way down to the faintest ISM emission. High 
latitude emission is all faint, and adds noise to the correlation at low intensities, but it still appears to follow 
the same correlation. 


6. Discussion 

TBD: The discussion will show how the constraints imposed by the new models can be used to estimate 
the 4.9 /xm CIB based on prior estimates of the 2.2 and 3.5 /xm CIB. Figure 8 shows how f all bands are 
constrained by various models and measurements, assuming a blackbody CIB function as an example. A 
more complicated CIB spectrum may be similarly depicted. Cautionary statements will be made about the 
similarity of the apparent CIB spectrum to the stellar and particularly the zodiacal light spectrum. 


7. Summary 

TBD: A summary, similar to the abstract, will appear here. Emphasis on the necessity for definite 
improvements in the zodiacal light modeling as the remaining major step in improving the detection of the 
near-IR CIB. 

This work was supported bv the NASA ADP program, grant NASW-99228. 


REFERENCES 

Arendt. R. G., et al, 1994, ApJ, 425, L85 
Arendt. R. G., et al. 1998, ApJ, 508, 74 
Dwek, E. & Arendt, R. G. 1998, ApJ, 508, L9 
Hauser, M. G., et al. 1998, ApJ, 508. 25 
Kelsall, T., et al. 1998, ApJ, 508. 44 



Rieke. G. H. & Lebofskv, M. J. 1985. Ap.J, 288. G18 


This preprint was prepareci with the AAS LAT^X macros v5.0. 



Fig. 1. — Images of 3.5 fi m residual emission after subtraction of foreground emission models, (a) Hauser 
et al. (1998), (b) Dwek k Arendt (1998), (c) Model 1, (d) Model 2. The images are shown as the Galactic 
azimuthal equal-area projections used by Hauser et al. with the north Galactic hemisphere on the left and 
south on the right. The Galactic center is at the bottom edge of each hemispherical image. The intensity 
ranges are linear from -0.01 to 0.05 MJv sr -1 for (a) and (b), and -0.0256 to 0.0344 MJy sr for (c). The 
bright source blanking that was applied to the Hauser et al. data was also applied to the other tw r o images. 
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Fig. 1. — Continued. 
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Fig. 2. — Histograms of the residual 3.5 jmi intensities in HQ AN and HQAS regions from (a,b) the Hauser 
et al. (1998) results, (c,d) the Dwek & Arendt (1998) results, (e,f) Model 1, and (g, h) Model 2. The smooth 
lines in panels (c) - (h) show Gaussian fits to the histograms. The histograms in (a) and (b) are fit with 
an additional parabolic base level, which mitigates the effect of the positive tail on the Gaussian fit to the 
peak. The HQ AN results are on the left, and HQAS results are on the right. 
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Fig. 3. — Same a 5 Figure 2. but for the smaller HQBN and HQBS regions. 












Fig. 4. — Gradients of the residual 3.5 /im emiss 
esc b. Panels (a,b) show the Hauser et al. (1998) i 
(g,h) show the results of Model 1 and Model 2. (i 
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Pig 5 — Two-point correlation functions of the 3.5 /ira residual emission calculated in the HQB regions for 
(a) the Hauser et al. (1998) result, (b) the Dwek & Arendt (1998) result, (c) the Model 1 result, and (d) 
the Model 2 result. The solid lines in each panel are the estimated ±la uncertainties. 



Fig. 6. — The top panel shows the 100 fxm ISM emission on a logarithmic intensity scale and a Mollweide 
projection. The lower panel shows the 3.5 ^ra ISM emission on a corresponding logarithmic scale, after 
subtraction of the Model 2 stellar emission. No bright source blanking has been applied. 
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Fig. 7. — ISM emission at 100 and 3.5 fim exhibits a tight linear correlation across a wide range of intensities. 
The intensities plotted here are from pixels at \b\ < 20°; higher latitude pixels add more data at the low 
intensity end of the plot, but do not alter the trend shown in the inset. The slope of the line drawn in the 
main panel and inset is the parameter C from Model 2 (Table 1). 
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Fig. 8. — Assuming a blackbody spectrum for the CIB, these panels show the various 2 a upper limits on the 
intensity at each wavelength, as a function of the CIB temperature. Upper limits from Hauser et al. (1998) 
are shown as dashed lines. Upper limits imposed by the constraints of Model 1 and Model 2 are shown by 
the black dashed and solid lines respectively. The Dw^ek k Arendt (1998) results imply both upper and lower 
limits on the CIB, which are shown as solid grey lines. 
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Table 1. Parameters of the 3.5 fin i Model 


Model Version 

A 

B 

c 

D 

Dwek k Arendt 

0.496 


0.00183 


Model 1 

0.445 

0.220 

0.00135 

-0.00106 

Model 2 

0.467 

0.250 

0.00122 

-0.00234 


Note. Values are for /„ in units of M.Iv sr 


Table 2. Systematic Uncertainties of the Mean 3.5 fin i Residual Emission 




Model 1 


Model 2 

Uncertainty 

a (nW m’ 2 

sr -1 ) Scale Factor 

a (nW in 2 

sr -1 ) Scale Factor 

Zodiacal Light (1.25 /xm) 



15 

0.0368 

Zodiacal Light (2.2 fim) 

6 

0.279 

6 

0.355 

Zodiacal Light (3.5 fim) 

2 

1.00 

2 

1.00 

Zodiacal Light (4.9 fim) 

6 

0.308 

6 

0.363 

Zodiacal Light (100 fim) 

6 

0.0387 

6 

0.0347 

Stellar &; ISM Model 

2.2 

1.00 

1.8 

1.00 

Total 

3.9 


4.1 



Table 3. Mean 3.5 fim Residual Emission at Select Patches 


Patch 

Size/Location 

Hauser et al. (1998) 

Dwek & Arendt (1998) 

Model 1 

Model 2 

NEP 

10° x 10° at d = +90 

5-6 

5.8 

-1.5 

— 1.1 

SEP 

10° x 10° at 3 - -90 

-2.3 

6.6 

-1.0 

-0.8 

NC.P 

10° x 10° at b = +90 

15.2 

11.4 

-0.1 

-1.6 

SOP 

10° x 10° at b = -90 

15.9 

12.7 

-0.3 

-2.1 

LH 

5° X 5° at (i, 6) = (150°, +53°) 

16.1 

11.0 

0.7 

-0.3 

HQ A 

\b\ > 30 and |0| > 25 

10.5 

9.9 

-0.1 

-0.9 

HQ AN 

b > +30 and 3 > +25 

11.3 

9.5 

-0.1 

-0.8 

HQ AS 

i6| < -30 and 3 < -25 

9.7 

10.2 

-0 0 

-0.9 

HQB 

[6] > 60 and |£i| > 45 

11.4 

9.9 

0.4 

-0.5 

HQBN 

b > +60 and 3 > -M5 

11.7 

9.5 

0.3 

-0.5 

HQBS 

£» < —60 and 3 < -45 

11. 0 

10.2 

0.4 

-0.5 


Note. — Uesults are for vl u in units of nW m 2 sr 
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Table 4. Gaussian Fits to 3.5 f.im Residual Emission Histograms 


Model Version 

Region 

Mean vl y 

(nW m~ J sr~ ] ) 

<y{vh) 

(nW m“ 2 sr ' 1 ) 

K-S Probability 
Non-Gaussian 

Hauser et al. (1998) 

HQAN 

5.70 

6.89 

1.00 

Hauser et al. (1998) 

HQAS 

4.25 

8.50 

1.00 

Dwek k Arendt (1998) 

HQAN 

9.35 

2.23 

1.00 

Dwek k Arendt (1998) 

HQAS 

10.09 

2.43 

1.00 

Model 1 

HQAN 

-0.09 

1.77 

1.00 

Model 1 

HQAS 

-0.05 

1.91 

1.00 

Model 2 

HQAN 

-0.78 

1.85 

1.00 

Model 2 

HQAS 

-0.90 

2.11 

1.00 

Hauser et al. (1998) 

HQBN 

5.05 

3.93 

1.00 

Hauser et al. (1998) 

HQBS 

3.96 

5.47 

1.00 

Dwek & Arendt (1998) 

HQBN 

9.44 

1.48 

0.93 

Dwek &; Arendt (1998) 

HQBS 

10.19 

1.74 

0.92 

Model 1 

HQBN 

0.31 

1.45 

0.83 

Model 1 

HQBS 

0.40 

1.63 

0.98 

Model 2 

HQBN 

-0.47 

1.50 

0.51 

Model 2 

HQBS 

-0.57 

1.71 

0.96 


Table 5. Galactic Gradients of the 3.5 /im Residual Emission 


Model Version 

Region 

Gradient 

(nW m~ 2 sr'^/cscdbl) 

Correlation 

Coefficient 

Hauser et al. (1998) 

HQA 

-9.91 ±0.23 

-0.16 

Dwek k Arendt (1998) 

HQA 

-3.32 ± 0.03 

-0.34 

Model 1 

HQA 

-0.96 ± 0.03 

-0.12 

Model 2 

HQA 

0.39 ± 0.03 

0.05 

Hauser et al. (1998) 

HQB 

-19.45 ±5.32 

-0.04 

Dwek k Arendt (1998) 

HQB 

-8.61 ±0.70 

-0.14 

Model 1 

HQB 

-2.53 ±0.66 

-0.05 

Model 2 

HQB 

1.39 ±0.67 

0.02 
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